The Hamiltonian limit of (3+l)D SU(3) lattice gauge theory on anisotropic lattices 
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The extreme anisotropic limit of Euclidean SU(3) lattice gauge theory is examined to extract 
the Hamiltonian limit, using standard path integral Monte Carlo (PIMC) methods. We examine 
the mean plaquette and string tension and compare them to results obtained within the Hamilto- 
nian framework of Kogut and Susskind. The results are a significant improvement upon previous 
Hamiltonian estimates, despite the extrapolation procedure necessary to extract observables. We 
conclude that the PIMC method is a reliable method of obtaining results for the Hamiltonian version 
of the theory. Our results also clearly demonstrate the universality between the Hamiltonian and 
Euclidean formulations of lattice gauge theory. It is particularly important to take into account the 
renormalization of both the anisotropy, and the Euclidean coupling Pe, in obtaining these results. 
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I. INTRODUCTION 

It is well known that lattice gauge theory (LGT) can 
be constructed in two ways: Wilson's "Euclidean" 
formulation where both space and time are discretized, 
or Kogut and Susskind's "Hamiltonian" Q formulation 
where time remains continuous. The method of choice 
in recent years for most lattice gauge theorists has been 
to use classical Monte Carlo methods in the Euclidean 
framework to extract observables from the theory. LGT 
in the Hamiltonian formulation has been rather neglected 
in comparison. Despite this, the Hamiltonian frame- 
work still offers an interesting alternative to its Euclidean 
cousin. One advantage of the Hamiltonian version is 
that techniques familiar from quantum many-body the- 
ory can be used to attack the problem, such as strong- 
couplingseries methods 's'l, the coupled cluster method 
B. IE the <-cxpansion 9, 10,, 11], the plaquette 

expansion [l2 . IT^ Hdj l , and the density matrix renormal- 
ization group (DMRG) method It must be said 
however, these methods tend to be more successful for 
lower dimensional lattices. Another possible advantage 
is that from a numerical point of view, the reduction in 
the dimensionality of the lattice from four to three pro- 
vides a significant reduction in computational overheads. 
Further, Hamiltonian results can serve as a check of the 
universality of Euclidean results. 

Due to the success of Monte Carlo methods in the Eu- 
clidean regime, one might expect similar levels of success 
in Hamiltonian LGT. Unfortunately, this has not been 
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the case and quantum Monte Carlo methods for Hamil- 
tonian LGT lag at least ten years behind their Euclidean 
counterparts. One of the first attempts at such a cal- 
culation was performed using a Green's function Monte 
Carlo approach by Heys and Stump for U(l) 0,0] and 
SU(2) JJi]. Chin and co-workers soon after implemented 
the closely relate d "g uided random walk" algorithm for 
U(l) in (3+l)D [Hi, then for SU(2) and SU(3) 

{21. 22. ^J. A feature of both these methods is that 
a "trial wave function" is used to guide random walk- 
ers towards the exact wavefunction. This appeared in 
the early days to be an advantage of the technique as 
the physical features of the wavefunction may be put in 
by hand, while the unknown part may be found through 
the stochastic process. However, later investigations [24l | 
have shown that there is an unacceptable dependence of 
the observables on the parameters of the wavefunction, 
in the case of SU(3) theory in (3+l)D. Other methods, 
such as the "projector Monte Carlo" method '25', "26^ and 
the related "stochastic truncation" method |27j have a 
version of the "minus-sign problem" arising due to the 
strong coupling (electric field) representation used for the 
basis states. This is due to the necessary introduction 
of Clebsch-Gordan coefficients for non-Abelian theories, 
which potentially cause destructive interference between 
the transition amplitudes. In view of the lack of any clear 
success of these quantum Monte Carlo methods, we are 
therefore forced to pursue an alternative approach. 

In a previous study, standard Euclidean path inte- 
gral Monte Carlo (PIMC) methods were used to extract 
the Hamiltonian limit for the U(l) lattice gauge theory 
in (2+1) dimensions The basic idea is to mea- 

sure observables on increasingly anisotropic lattices, then 
extrapolate to the Hamiltonian limit, corresponding to 
At = at/ as — > 0, where at and are the lattice spacings 
in the time and space directions respectively. The results 
obtained show excellent agreement between the ex- 
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trapolated results and Hamiltonian estimates of the same 
quantities. In this work we attempt a similar procedure 
for pure SU(3) gauge theory in (3+1) dimensions. We 
calculate two basic quantities, the average plaquette and 
the string tension, and compare them to results obtained 
in the Hamiltonian formulation. We find that indeed the 
anisotropic Euclidean results converge to the Hamilto- 
nian estimates, once the difference of scales between 
the two theories has been taken into account. Specifi- 
cally we find that this is particularly important at finite 
anisotropy |30l | in the extrapolation procedure. This will 
be discussed in more detail in section HTl 

The use of anisotropic lattices has generated much 
interest in recent years, following the work of Morn- 
ingstar and Peardon |^ [s^ in extracting accurate glue- 
ball masses. This same approach has now been extended 
to extracting heavy quark spectra [ssL Is^ . The impor- 
tance of taking into account the difference of scales on 
anisotropic lattices in these studies has already been dis- 
cussed by Klassen [s^, through the renormalization of 
the anisotropy. In our case, the renormalization of the 
the isotropic Euclidean coupling f3E is also important in 
extracting our results. This may be of some relevance to 
other studies on anisotropic lattices. 

In Sec. ^we briefly discuss the SU(3) model as defined 
on anisotropic lattices. We also discuss in this section our 
extrapolation procedure to the Hamiltonian limit. In Sec. 
mil we explain the simulation methods used to obtain our 
results, and in Sec. II VI we show our results for the average 
plaquette and string tension. Some concluding remarks 
are given in Sec. El 



II. THE ANISOTROPIC SU(3) MODEL 

The anisotropic action is given by [3^ 

X i>j- i^A X i^A 

(1) 

where the first term sums over all space-like plaquettes 
on the lattice, and the second term sums over time-like 
plaquettes. A plaquette at lattice position x is defined 
by 

P,,,{x) = ^ReTr [U^,{x)U,{x + /i);7^(x + v)Ul{x)\ , 

(2) 

where U^{x) is the SU(3) gauge field variable. The cou- 
plings Pa and j3r are defined by 



-- 



B _1 
9t 



(3) 



and the anisotropy factor, or aspect ratio is defined by 

1 a. 



One must include different space-like and time-like cou- 
plings ga and gr in Eq. Q in order to allow the free- 
dom to renormalize so that correlation lengths are equal 
in both directions, even though the spacings as and at 
arc different. In the continuum limit we require that 
physical quantities be independent of changes in both of 
these quantities. Therefore we require the two couplings 
to be a function of both as and «(. Karsch ^30] , u sing 
the background field method of Dashen and Gross |33 
and Hasenfratz and Hasenfratz pof , obtained a mapping 
between the equivalent couplings of the Euclidean and 
anisotropic actions 
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5f 



9l 



\+c,{i) + 0{gl), 
9e 

^ + c.(e) + o(5i), 

9e 



(5) 
(6) 



where is the coupling for the Euclidean theory with 
^ = 1. The factors c„{^) and Cr{S,) are defined in Eqs. 
(2.24) and (2.25) of Ref. [13 . For the Euchdean theory 
with = 1; these factors both approach zero, and there- 
fore we recover the usual action with only one coupling 
.9ct = gr = gE- In the limit ^ —t 00, Eqs. Q and ® 
reduce to the Hamiltonian values as obtained by Hasen- 
fratz and Hasenfratz 29]. 

It is convenient to rewrite the action in a more 
symmetric fashion 



/5? 



where 



with 



^ J2 (i-p,,(x))-i-/3^eEE(i-^4.(2;)), 



X i>j; 



X i^A 



B - ^ 

Pi — —y 



(7) 



(8) 



5{ = 5<t5t and 77 = [gl/glf/'^. In the limit ^ ^ 00, 
goes to the Hamiltonian coupHng fin- Using Eqs. © 
and © we may write a correspondence between these 
variables and the Euclidean theory 



= /3£ + 3(c,(0 + c,(0) + O(/3£i 



Pe 



Cr(0) + O(/3/) 



(9) 
(10) 



At at 



(4) 



where Pe — ^l9%- Therefore for every {Pe,£,) pair there 
is a corresponding pair of couplings (/3^,0- The relation 
between £ and f has been discussed in some detail by 
Klassen [33. In his language, f is the bare anisotropy, 
while £ is the renormalized anisotropy. To evaluate the 
factors Cct(£) and Ct-(£), one may either directly calculate 
them in terms of the integrals given in Ref. 30], or in 
the case of Eq. (|10|l , use the parameterization given by 
Klassen [ssl ]. 

We now discuss our extrapolation procedure in order 
to obtain Hamiltonian estimates from an anisotropic lat- 
tice. In a naive extrapolation procedure, one might as- 
sume P = Pa- = Pr in Eq. (Q, and extrapolate phys- 
ical quantities at constant P to the Hamiltonian limit. 



0.2 




FIG. 1: The discrepancy between the anisotropic and Eu- 
chdean couphngs as a function of the anisotropy. 

^ — > oo. This procedure is incorrect, however, because 
/3ct 7^ /3r 7^ /3 due to renormalization. The correct pro- 
cedure is to extrapolate to ^ ^ oo at constant Pe- We 
may summarize our procedure as follows: 

1. For a particular /J^, choose several anisotropics ^. 

2. Calculate the corresponding values of and ^ us- 
ing Eqs. lO and ^U^. 

3. Use the action given in to calculate physical 
observables at these couplings. 

4. Perform a polynomial fit to the data in the inverse 
square anisotropy At^ at constant and extrap- 
olate measured quantities to ^ ^ cxd (At — > 0) (we 
extrapolate in the inverse square anisotropy since 
the leading discretization errors in the timelike di- 
rection are expected to be of order Oj .). 

It is interesting to note that the discrepancy between 
and Pe in Eq. Q reaches a maximum around f « 3.6215, 
as can be seen from Fig. Q] We see that in the vicinity of 
the maximum, the discrepancy reaches nearly P(_ — Pe ~ 
0.19. Since the most anisotropic lattice we use is in the 
vicinity of this maximum, significant discrepancies begin 
to creep in if the correct procedure is not used. 

III. METHOD 

A. Simulation Details 

We analyze the action given in Eq. ^ by stan- 
dard path integral Monte Carlo methods. Configurations 
are generated using the Cabibbo-Marinari pseudo- 
heat bath algorithm, applied to the three diagonal SU(2) 
subgroups of SU(3). To analyze the behavior between 
weak and strong coupling we concentrate on the region 
/3j = 4 — 7. For each P^ value we generate configurations 
for anisotropics in the range ^ = 1 — 3. These parameters 
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TABLE L Parameters used for each configuration set. 



Dimensions 


ft 




Pe 




8^ X 8 


4.0 


1 


4.0 


1.0 


8^ X 8 


5.0 


1 


5.0 


1.0 


8^ X 8 


5.4 


1 


5.4 


1.0 


8^ X 8 


5.6 


1 


5.6 


1.0 


8^ X 8 


5.8 


1 


5.8 


1.0 


8^ X 8 


6.0 


1 


6.0 


1.0 


8^ X 8 


7.0 


1 


7.0 


1.0 


8^ X 12 


4.0 




3.9396 


1.2968 


8^ X 12 


5.0 


^3/2 


4.9428 


1.2795 


8^ X 12 


5.4 


^3/2 


5.3437 


1.2747 


8=* X 16 


4.0 


V2 


3.8992 


1.5444 


8^ X 16 


5.0 


V2 


4.9036 


1.5139 


8^ X 16 


5.4 


V2 


5.3049 


1.5053 


8^ X 12 


5.0 


3/2 


4.8891 


1.6197 


8^ X 12 


5.4 


3/2 


5.2904 


1.6095 


8^ X 12 


5.6 


3/2 


5.4910 


1.6050 


8^ X 12 


5.8 


3/2 


5.6916 


1.6008 


8^ X 12 


6.0 


3/2 


5.8921 


1.5970 


8^ X 12 


7.0 


3/2 


6.8941 


1.5815 


8^ X 24 


4.0 


V3 


3.8546 


1.9569 


8^ X 24 


5.0 




4.8588 


1.9053 


8^ X 24 


5.4 


V3 


5.2602 


1.8907 


8^ X 16 


5.0 


2 


4.8366 


2.2341 


8^ X 16 


5.4 


2 


5.2376 


2.2146 


8^ X 16 


5.6 


2 


5.4381 


2.2060 


8^ X 16 


5.8 


2 


5.6385 


2.1981 


8^ X 16 


6.0 


2 


5.8389 


2.1908 


8^ X 16 


7.0 


2 


6.8406 


2.1609 


8^ X 24 


4.0 


3 


3.8110 


3.5811 


8^ X 24 


5.0 


3 


4.8112 


3.4538 


8^ X 24 


5.4 


3 


5.2113 


3.4171 


8^ X 24 


5.6 


3 


5.4114 


3.4010 


8^ X 24 


5.8 


3 


5.6114 


3.3860 


8^ X 24 


6.0 


3 


5.8115 


3.3721 


8^ X 24 


7.0 


3 


6.8117 


3.3152 



{P^, ^) may be converted into their corresponding {PetO 
by solving © and (|10|l numerically. Our full set of pa- 
rameters is shown in Tabled together with the equivalent 
{Pe.Ci values. 

Due to the large number of configuration sets that must 
be generated, we limit ourselves to relatively modest lat- 
tice sizes, ranging from 8'^ x 8 to 8'^ x 24. We adjust 
the lattice size in the time direction according to the 
anisotropy used in order to keep the physical length in 
the time and space directions equal. Configurations are 
given a cold start, and then 5000 thermalization sweeps 
in order to equilibrate. After thermalization, configura- 
tions are stored every 500 sweeps, for 100 configurations. 
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5 = 1.0 
^=1.5 
1=2.0^ 
^ = 3.0 ^ 



20000 30000 

N„„..„ 



were constant in Pe, by choosing appropriate varying 
(/^fiO values. In a future study we would adopt such a 
procedure, although there is little loss in accuracy for our 
current results. 



B. The Static Quark Potential 

The static quark potential is extracted from the ex- 
pectation values of Wilson loops, which are expected to 
behave like 



M/(r,<)=^a(r)exp(-l^,(r)i), 



(11) 



FIG. 2: Evolution of the average action (S) with the number 
of sweeps of the lattice, for f3^ — 5.6, and various anisotropics. 



As a measure of the equilibration, we plot the average 
action for various anisotropics in Fig. |21 We see that for 
each value of the anisotropy the configurations relax to 
equilibrium after of order 1000 sweeps. In particular we 
see very little difference in equilibration times for the var- 
ious anisotropics. This is in contrast to the results of Ref. 
[2^, which used a standard Metropolis algorithm com- 
bined with Fourier acceleration techniques for the U(l) 
model in (2+1) dimensions. There it was found that 
equilibration times were longer for the anisotropic lat- 
tices, despite using Fourier acceleration techniques. The 
advantage of the pseudo-heat bath algorithm in updating 
anisotropic lattices, where the space-time (3t and can 
be very large, is most easily seen when considering its ap- 
plication to SU(2)-color gauge theory. In this case a link 
variable is updated without any reference to the current 
link. Rather only the staples associated with the link are 
considered in making the SU(2) update. The improved 
ergodicity afforded by this algorithm is largely carried 
over to the SU(2) subgroup updates of SU(3) links. 

According to the extrapolation procedure described 
in Sec. ^ we must extrapolate at constant Pe, not 
/3j as we have calculated. Therefore for all observables 
we must interpolate our values in the (/3^,0 plane onto 
lines of constant Be- We do this using the modified 
Shepard method [39j of multivariate interpolation. We 
choose our interpolation points along the lines (3e — 
3.8110, 4.8112, 5.2113, 5.4114, 5.6114, 5.8115, 6.8117, and 
^ values such that they are as close as possible to existing 
data points. These points are chosen such that there is no 
interpolation for the most anisotropic points = 3), and 
there is minimal interpolation for the remaining points, 
to minimize the interpolation error. In practice these in- 
terpolation errors are quite negligible compared with the 
original errors on the data points, such that the error on 
the interpolated value can be taken as the same as the 
error in the nearest data point. 

Clearly a more straightforward procedure would have 
been to generate data points originally such that they 



where the summation is over the excited state contribu- 
tions to the expectation value, and i = I corresponds to 
the contribution from the ground state. To obtain the 
optimal signal-to-noise ratio, we must suppress the con- 
tributions from the excited states, which may be done by 
APE smearing [iol l4lj | . This involves replacing a partic- 
ular space-like link by 



+Ul{x - £')Uf,{x ^ i>)U^{x ~ v + fl)] , 

(12) 

where the P denotes a projection back onto SU(3), and 
a is a parameter which is defined by the user. This is 
repeated n^pE times. This amounts to adding in a frac- 
tion a of neighboring staples for every space-like link on 
the lattice. To tune the smearing parameters, we fix the 
smearing fraction to a = 0.7, as it is sufficient to fix a 
and tune riAPE To find the optimum tt-apEj we look 
for the cleanest plateaus in the effective potential 



V{r, t) = In 



_W{r,t + l) 
and also examine the ratio [43l | 

T4^*+i(r, t)/W\r,t + l), 



(13) 



(14) 



which should be near unity for good ground state domi- 
nance. A typical value which proved to be sufficient for 
most cases was ^ape — 5, although for small the over- 
all effectiveness of the smearing procedure was reduced. 
An example of a typical effective potential plot is shown 
in Fig. 13 for (3^ = 5.60 and f = 2.0. We see good plateau 
behavior for small t values, which reflects the optimum 
smearing. Finally we extract the string tension by fitting 
with the form 



V{r) = Vo + Kr-e/r, 



(15) 
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FIG. 3: Effective potential plots for 13^ = 5.60 and ^ = 2.0. 



V(r) 




1 2 3 4 5 6 

r 



FIG. 4: Static quark potential at f3^ = 5.8 and ^ = 3.0. 

where e = 7r/12 0], and K and Vq are fit parameters. 
The string tension is then found according to K = era?. 
An example of such a fit is shown in Fig. 01 for /3j = 5.8 
and ^ = 3.0. We see that the data are fitted very well 
by (|15|) . in this instance giving K = aa} = 0.189(4) and 
Vo = 0.732(7). 



IV. RESULTS 

A. Average Plaquette 

Our aim in this section is to obtain estimates for the 
average plaquette, defined as the expectation value of Eq. 

for the Hamiltonian theory. As no time-like plaque- 
ttes are present in the Hamiltonian lattice, we therefore 
only calculate space-like plaquettes in the Euclidean the- 
ory. Henceforth we will refer to the average space-like 
plaquette simply as "average plaquette" , unless specified 
otherwise. As a first test of our results, we compare the 
average plaquette for the isotropic lattice ^ = 1 to exist- 



FIG. 5: Comparison of our PIMC data for the average 
plaquette with strong and weak coupling expansions for the 
isotropic lattice ^ = 1. Dashed lines are merely to guide the 
eye. 

ing results. This is shown in Fig. | 5 l to gether with strong 
coupling expansions to order /3 1? [4^ |4a| , and weak cou- 
pling expansions to order jig. We also show a [5/5] 
Fade extrapolation of the strong coupling series, which 
shows that beyond (3e ~ 5 the series diverges away. Our 
PIMC data match smoothly onto the strong coupling and 
weak coupling expansions in their respective limits, as ex- 
pected |58| . It can be seen that the crossover from strong 
coupling to weak coupling behavior takes place in the re- 
gion « 5 — 6, as is well known. 

At finite anisotropy, we are not aware of any results 
available to make a similar comparison. It is fairly 
straightforward however to calculate the corresponding 
strong coupling expansion to order The calculations 
for this are given in the Appendix. Setting I = J = 1 
into (|A.10|I gives 

where b = /3^/6. A comparison between the strong cou- 
pling expansion and our PIMC data is shown in Fig. El 
along with a [2 /2] Fade extrapolation. We see consistent 
results between the expansions and our data for f3^ = 4, 
which is still in the strong coupling region. There is a 
trend for the series with larger values of the anisotropy 
^ to break down earlier as can be seen from the Fade 
approximant for ^ — y/3. 

We now extrapolate our results to the Hamiltonian 
limit. Performing the interpolation as described in Sec. 
nil Al we obtain points at constant Pe- Our results are 
then extrapolated to the Hamiltonian limit in powers of 
Ar^, as shown in Fig. [7| We see a fairly smooth de- 
pendence on Ar^, for all (3e- Error estimates for the 
extrapolation may be estimated by comparing fits of dif- 
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TABLE II; A summary of our estimates in the Hamiltonian 
limit Ar — > for the average (space-hke) plaquette and the 
string tension. The Hamiltonian coupling A, calculated from 
Eq. 1171 . is also shown. 



Pe 


A 






3.8110 


2.5214 


0.125(5) 




4.8112 


3.9849 


0.23(1) 


1.08(20) 


5.2113 


4.6637 


0.312(7) 


0.57(4) 


5.4114 


5.0231 


0.36(1) 


0.31(2) 


5.6114 


5.3958 


0.403(5) 


0.17(1) 


5.8115 


5.7819 


0.431(6) 


0.110(4) 


6.8117 


7.9125 


0.526(7) 





FIG. 6: Comparison of the average space-like plaquette cal- 
culated using PIMC at anisotropies f — -^3/2 and — \/3 
with strong coupling expansions to order /3|. Dashed lines are 
to guide the eye. 
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P^4.81 
Pj=5.21 



Pj=5.41 ^•^ 
P^5.61 ^B- 
Pj=5.81 
Pi,= 6.81 ^X- 



series to order A^^/^ of Hofsass and Horsley ^'^d 
strong coupling series to order A*" of Hamer, Irving and 
Preece For the strong coupling series we show both 
the raw series sum, and the [3/3] Pade approximant. The 
analysis of Hamer et al. 0| shows that finite-size effects 
in the Monte Carlo data should be negligible at this level 
of accuracy. 

For A < 4, the agreement between the different esti- 
mates is excellent, except that our point at A = 2.52 is 
a little lower than the strong coupling series result. A 
possible explanation is that we have relied heavily on the 
weak-coupling relations between scales of the anisotropic 
and Hamiltonian theories and the Euclidean theory (i.e. 
Eqs. (O, ((Tn|) . and ifTTjl ). It is possible higher order terms 
in these expansions begin to contribute larger effects at 
A < 3. 



FIG. 7: Extrapolation of the average plaquette to the Hamil- 
tonian limit Ar — > 0, for fixed I3e- Dashed lines show 
quadratic fits to the data in Ar^. 

ferent orders (linear, quadratic, cubic, ...) in Ar^. Our 
numerical values in the Ar — > limit are shown in Table 

m 

To compare our extrapolated results to existing Hamil- 
tonian results, we must take into account the difference 
of scales between the two regimes. The Hamiltonian cou- 
pling parameter A — 6/gfj, where gn = lim^^oo^j, may 
be related to the Euclidean coupling through the relation 

M 

= \/6A- 0.07848. (17) 

Using this relation wc may plot our extrapolated 
Euclidean estimates against previous results obtained 
within the Hamiltonian formulation, as shown in Fig. 
|S1 This figure compares our present results with Green's 
Function Monte Carlo (GFMC) results of Chin et al. [U 
and Hamer, Samaras, and Bursill j2J|, the weak coupling 



In the region A = 5 — 8, however, there is a large dis- 
crepancy between our results and the previous GFA^IC 
estimates jj^, i2M ■ We find that the crossover from strong 
to weak coupling behavior occurs at much the same cou- 
plings as in the isotropic regime, and earlier than shown 
by the GFMC results. There is quite good agreement, 
on the other hand between our results and the strong 
coupling series extrapolations in this regime. This gives 
us confidence that our present results are the more ac- 
curate. It was pointed out in Ref. (231 that the GFMC 
method suffers from an unacceptable dependence on the 
trial wave function, and we surmise that the variational 
parameter may not have been optimized to the correct 
value in this region. 

Further evidence is provided by Long et al. (23| who 
showed that GFMC results for the ground-state energy 
change significantly if a second variational parameter is 
added in to the trial wave function. We conclude that 
the PIMC offers a more unbiased method of extracting 
results, even with the somewhat undesirable extra step 
of having to perform a At — > extrapolation. 
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Chinetai. 
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[3/3] Fade to strong coupling 




FIG. 8: Our extrapolated Euclidean PIMC results compared 
to Hamiltonian strong coupling 49], Hamiltonian weak cou- 
pling 13, and Green's Function Monte Carlo (GFMC) results 
|24|. Dashed lines are merely to guide the eye. 



B. String Tension 

We now turn to calculating the string tension. We 
again use only space-like Wilson loops, so we can ulti- 
mately compare to Hamiltonian results. Fig. O shows 
results for the isotropic case, ^ = 1. Our results for small 
(3e are contaminated by large errors due to the fast ex- 
ponential decay of the Wilson loops. The parallel lines 
show the two-loop scaling form 



1 



where 



70 



11 

16^ 



71 



102 

(167r2)2' 



(18) 



(19) 



The coupling g refers to the Euclidean coupling gs in 
this case, and the scaling parameter here is A = A^;, the 
Euclidean scaling parameter. Our data asymptotically 
appear to approach the expected scaling form in the weak 
coupling region. The parallel lines correspond to 



(113± 10)A£;, 



y/a^ (1.35±0.12)A mom: 



(20) 



(21) 



where Amom is the perturbative QCD scale parameter 
|29l l50j |. It is well known, however, and is evident from 
Fig. 1^ that in quantitative terms the string tension fol- 
lows neither two-loop or three-loop scaling at these cou- 
plings; and a more sophisticated fit by Edwards et. al. 
|5l| taking into account various correction terms gives a 
much smaller value. 



(22) 



FIG. 9: The string tension at ^ = 1. We plot our own PIMC 
results, together with a 12th order strong coupling expansion 



We do not attempt such a fit here. 

On the same figure are plotted the strong cou plin g ex- 
pansions of Miinster and Weisz to 12th order [SjI, and 
the Monte Carlo results of Edwards et al. [HJl . Due to 
the presence of the roughening transition at (3e ~ 5.5, we 
expect the strong coupling expansion to diverge from the 
data beyond this point, as is seen in the figure. At the 
smaller values of /3_Bj_our results agree well with those 
of Edwards et al. but at larger (3e, our results 

are higher than theirs. This is most probably due to the 
smaller lattice size used in our calculations, as discussed 
in Ref "sTl , and is consistent with the higher estimate of 
the string tension obtained at equation H21|l . One would 
need to use larger lattice sizes or else an improved ac- 
tion to get more accurate estimates of the asymptotic 
parameters. 

The extrapolation to the Hamiltonian limit is shown 
in Fig. 1101 This is performed again in powers of Ar^, 
in a similar fashion to the average plaquette. There is a 
trend towards a stronger curvature in the extrapolation 
for smaller values of (5ei suggesting that our estimate for 
I3e = 4.81 is probably too high. There are however large 
error bars on these data points. 

In Fig. our extrapolated results for the string ten- 
sion are compared with Hamiltonian estimates obtained 
by Hamer, Irving and Preece using an Exact Linked 
Cluster Expansion (ELCE) method. Our axes are such 
that the scaling relation (|18|) for the string tension ap- 
pears as a straight line in the plot. In this case we use the 
Hamiltonian versions of the couplings, g = gn, A = Ag, 
and A = G/gj^. We do not show here the GFMC results 
of Hamer et al. '2^, which were derived from Creutz 
ratios on very small loops, and therefore subject to large 
finite-size effects. We see that the ELCE results are more 
accurate in the strong coupling region, but our PIMC re- 
sults are better in the weak coupling region. The two 
sets of data do not agree below VA = 2.1. This may be 
because our point at VX = 2.00 (or (3e = 4.81) is too 
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FIG. 10: Extrapolation of the string tension to the Hamil- 
tonian limit At — > 0. Dashed lines show quadratic fits in 



much more reliable than the GFMC results of Ref. [2J| , 
mentioned above. 

We again see evidence for the approach to asymptotic 
scaling behavior in the region \f\ > 2.3. The data give 
us a fit of 



Vct = (123±5)A//. 



(1.34±0.05)Amom, 



(23) 



(24) 



using the conversion factors computed by Hasenfratz and 
Hasenfratz This is much the same as our Euclidean 
estimate, and somewhat too high by comparison with Eq. 
(|22|l . presumably for the same reasons discussed above. 
It is considerably better than the Hamiltonian ELCE es- 
timate ^/a = (1.7 ± 0.4)Aniom of Hamer et al. 49], how- 
ever. 




FIG. 11: Extrapolated PIMC results together with ELCE 
results, and GFMC estimates of the string tension. 



high, as noted above; but also, we do not expect precise 
agreement in this region, since the measurements corre- 
spond to different quantities. The PIMC estimates refer 
to the decay exponent of space- like Wilson loops, while 
the ELCE estimates are of the 'axial' string tension, i.e. 
the energy per unit length of a 'string' of flux along one 
axis. It is well known that these estimates may differ at 
strong coupling, where rotational invariance is broken; 
but at some intermediate coupling estimated [ 5 3 to be 
around A = 5.4, a "roughening" transition |44l l54 Issj 
takes place, after which rotational invariance is restored 
and different estimates of the string tension should coin- 
cide. Series expansions for the axial string tension cannot 
be continued past the roughening transition. 

Again the crossover to weak coupling behavior occurs 
in the region ^/X — 2.1 — 2.4, corresponding to (3e — 
5.1 — 5.8. In the weak coupling region, the ELCE results 
are somewhat higher than our PIMC estimates, but are 
almost compatible within errors: we believe our present 
results are more reliable and accurate. They are also 



V. CONCLUSIONS 

We have demonstrated in this work that one can obtain 
reliable results for the Hamiltonian limit using the stan- 
dard Path Integral Monte Carlo method for anisotropic 
lattices and extrapolating to At 0. We have cal- 
culated two quantities, the average (space-like) plaque- 
tte and the string tension. In both cases our results 
are a substantial improvement on other estimates cal- 
culated purely within the Hamiltonian formulation, al- 
though there is broad agreement between them. This 
demonstrates clear evidence of the universality between 
the Hamiltonian and Euclidean formulations of LGT. 

In performing the extrapolation to the Hamiltonian 
limit we have found it very important to take into ac- 
count the renormalization of the couplings /3j and ^ in 
Eq. ITJ at finite ^. The renormalization of £ has been 
discussed in some detail already by Klassen [23 ■ In our 
case this is only half the story, as there is also the renor- 
malization of (3(_, which is related to the Euclidean cou- 
pling /5b, through (Q. As mentioned in the text, the 
discrepancy between the two couplings reaches a maxi- 
mum around £ = 3.6214, which is quite close to some 
of our most anisotropic data points, and has a large in- 
fluence in terms of extrapolating our results to £ ^ oo. 
We expect that this renormalization will have a much 
smaller effect for improved actions [5^ [H^], and hence 
utilizing such actions may be a neater way of removing 
the complication altogether. 

Our main motivation for examining anisotropic lattices 
was to investigate alternative Monte Carlo procedures for 
obtaining results in the Hamiltonian formulation, in view 
of the lack of a reliable quantum Monte Carlo algorithm 
for Hamiltonian LGT. We have found that the PIMC ap- 
proach is superior to quantum Monte Carlo methods, in 
particular the GFMC algorithm. The GFMC algorithm 
was found previously |24i | to have unacceptable system- 
atic errors due to its dependence on a trial wave func- 
tion. In this investigation we have found that the PIMC 
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method gave more reliable results for the mean plaque- 
tte, and also gave good results for the string tension in 
the scaling regime, which the GFMC method |23| could 
not. 

The major disadvantage of PIMC is the necessity to 
make an extrapolation to the Hamiltonian limit Ar 0. 
This reduces the accuracy of the final results, and also is 
rather expensive in computer time, as several configura- 
tion sets must be generated to obtain one (3e point. With 
modern algorithms, however, the results are almost as ac- 
curate at large anisotropics as in this isotropic case, and 
the extrapolation can be made with confidence - see Figs. 
ITIandlTUI Thus we conclude that PIMC is the preferred 
Monte Carlo approach for estimating physical quantities 
in the Hamiltonian limit, just as in the Euclidean case. 

An obvious extension of the present work would be to 
attempt to estimate glueball masses using the PIMC ap- 
proach, and to compare those with existing Hamiltonian 
results. We hope to attempt this in future work. 



The coefficients dr and fr may be found by inverting the 
relation 

dr = J[dU]x:{U)e-'''^''\ (A.6) 

and similarly for fr- We only need the result for r = 
and r = 3 however, which may be found in Eqs. (7) and 
(10) of Ref. 0. Defining 



u{x) 



1 2 ^4 ^ 
X + —X — —X — —X 



7 



(A.7) 



72 24 720 
then coefficients 63 = d^/do and = /3//0 are given by 



53 



3u(/3jf76) 



6^ ■ 72^2 
6 72 



(A.8) 
(A.9) 



APPENDIX: STRONG COUPLING EXPANSION 
FOR SU(3) ON ANISOTROPIC LATTICES 

A Wilson loop of size / x J is defined by 



W{C) 



\ I m^x ( n u.) 



n 



X I I e 

PeP, 



-S"{P) 



n ^ 

PePst 



-S=*(P) 



(A.l) 



where C denotes the contour of the Wilson loop, N is 
the dimension of the group matrices (in our case A^ = 3), 
Uij denotes a group element lying between sites i and 
j, Z is the normalizing factor, [dU] denotes the group 
integration over SU(3) matrices, xO takes character of its 
argument, the product over P^s runs through all space- 
like plaquettes, while Pgt denotes all time-like plaquettes. 
We also define 



xiu) 



(A.2) 
(A.3) 



Using standard character expansion techniques [Sg, we 
write 



^drXr{U) 
frXriU). 



00 

E 

r=0 



(A.4) 
(A.5) 



Performing expansions according to the diagrams de- 
scribed in Ref. ;363 , we obtain for space- like Wilson loops 



~2IJN I 



21JN { — 



(A.IO) 



while for time-like Wilson loops we obtain 



(I 



.2 /ggN 3 



.N 



(A.ll) 



The string tension may then be computed using these 
results via aa^ ^ (1//./) \nW{C). 
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tices is sufficiently small to take us beyond the deconfine- 
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extent constrains us to examine short distance quantites 
which are not severely affected by deconfinement issues. 
Therefore, we present results for the mean plaquette at 
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